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ABSTRACT 


The response of an insulated wall, over which a heated fluid flows, to a sinusoidally 
forced fluid temperature was used to calculate the convective heat-transfer coefficients. 
An exact solution is given which accounts for thermal conductivity and the location of the 
sensed wall temperature in one-dimensional heat-transfer problems. Charts are in- 
cluded to aid in the calculation. A comparative analysis was made of solutions that do not 
account for thermal conductivity and the location of the sensed wall temperature and those 
that do. If the exact solution is not used, errors greater than 23 percent are possible. 
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DETERMINATION OF CONVECTIVE HEAT-TRANSFER COEFFICIENTS 
ON ADIABATIC WALLS USING A SINUSOIDALLY FORCED 
FLUID TEMPERATURE 
by Ronald G. Huff 
Lewis Research Center 

SUMMARY 

The response of an insulated wall, over which a heated fluid flows, to a sinusoidally 
forced fluid temperature was used to calculate the convective heat-transfer coefficients. 
An exact solution is given which accounts for thermal conductivity and the location of the 
sensed wall temperature in one-dimensional heat-transfer problems. Charts are in- 
cluded to aid in the calculation. A comparative analysis was made of solutions that do not 
account for thermal conductivity and the location of the sensed wall temperature and those 
that do. If the exact solution is not used, errors greater than 23 percent are possible. 


INTRODUCTION 

Both transient and steady-state analyses have been used to design calorimeters for 
use in rocket engines and aerodynamic heat-transfer studies. The steady-state calorim- 
eter makes use of the temperature gradient in a material of known conductivity and geom- 
etry; the transient type makes use of the response of a material to a driving temperature, 
that is, the response of a thin disk to a step change in the surrounding fluid. 

The response of a wall to a fluid that flows over it and has a sinusoidally oscillating 
temperature has been used to calculate the convective heat-transfer coefficient h. The 
solution found by Anderson (ref. 1) for the convective heat -transfer coefficient (herein 
called slug solution), however, does not account for thermal conductivity or the location 
of the measured wall temperature. An estimate, given by Anderson, of "the error in the 
measured time constant r caused by heat conduction through the skin” is 1. 5 percent. 
Bell (ref. 2) neglects the effect of thermal conductivity by designing his experiments in 
such a way as to cause its effect to drop out of his equations, which are in series form. 



The objective of this analytical investigation, conducted at NASA Lewis Research 
Center, was to find a solution for the convective heat-transfer coefficient as a function of 
the phase lag between the fluid and wall temperatures. Such a solution would take into 
account the thermoconductivity as well as the location of the measured wall temperature. 

The solution is presented along with charts that (for the wall temperature measured 
at the insulated side of the wall) can be used to determine the heat-transfer coefficient as 
a function of frequency, wall properties, wall thickness, and phase lag. A comparison is 
made between this solution and that of Anderson (ref. 1), both of which assume that the 
back surface of the wall x = L is perfectly insulated. 

The slug solution may be substituted for the present analysis when conductivity and 
temperature-sensor location are not important. 


SYMBOLS 

C specific heat of wall material, Btu/(lb)(°R); J/(kg)(K) 

CON function defined in eq. (4b) 

f frequency of temperature oscillation, cps; Hz 

h convective heat-transfer coefficient on surface of wall, Btu/(in. ^)(sec)(°R) ; 

W/(m 2 ) (K) 

K thermal conductivity, Btu/(in.)(sec)(°R); J/(m)(sec)(K) 

L thickness of wall, in. ; m 

T wall temperature, T w - Tq, °R; K 
T mean temperature, °R; K 

ATq amplitude of gas temperature, Tq - Tq, °R; K 
x distance measured from fluid-wall interface into wall, in . ; m 

or thermal diffusivity, K/pC 

£q function defined after eq. (4f) 

7] frequency and diffusivity perimeter, Vcu/ 20 1 

9 time, sec 

n constant equal to 3. 1416 rad 

O O 

p density of wall material, lb/in. ; kg/m 

r time constant, pCL/h 
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<p phase shift between forced fluid temperature Tq and wall temperature T , < 0 for 
T w lagging T Q , deg 

co angular velocity of forced fluid temperature Tq, rad/sec 
Subscripts: 

G fluid flowing over wall 

s values calculated with Anderson's slug-type solution (ref. 1) 
w wall over which fluid flows 

DIFFERENTIAL EQUATIONS AND ASSUMPTIONS 

The solution for the heat-transfer coefficient as a function of phase lag <p, frequency 
of fluid-temperature oscillation f, wall properties, and location of the sensed wall tem- 
perature T is now given. Consider an infinite plate on one side of which a fluid flows 
over the surface (x = 0) . The other side is insulated (x = L) . The system is illustrated 
in figure 1 . 


Fluid flow 


Surface exposed to fluid __ 

y//////A^m surfac rZZZZW, “ 
Figure 1. - Basic heat-transfer model. 


The applicable differential equation is 


o 2 T _ 1 3T 

ax 2 " “ 39 


where 

T difference between wall temperature at any location in wall and average fluid tem- 
perature 

x distance into wall from fluid side of wall 
a thermal diffusivity, K/pC 
0 time, sec 
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This solution assumes that the 

(1) Thermal conductivity is finite and constant 

(2) Convective heat -transfer coefficient h is constant 

(3) Density p and specific heat C of the wall are constant 

(4) Surface of the wall at x = L is insulated 

(5) Surface of the wall at x = 0 is exposed to a fluid whose temperature is given by 

T q - T q = AT G e' ia,e (2) 

where 

T g fluid temperature 

T g average fluid temperature 

AT g amplitude of gas temperature 

co angular velocity of temperature oscillation, 2?rf rad/sec 

(6) Convective heat transfer at x = 0 is 

-K — = h r (T r - T) (3) 

0x G G 

(7) Heat conduction through the wall is one dimensional 

For the solution to the differential equation (1), a product solution is assumed and the 
boundary conditions are applied (i.e. , an insulated surface at x = L, assumption (4), and 
convective heat transfer at x = 0, assumption (6) . The details of the solution are given 
in the appendix. The convective heat-transfer coefficient is 

hp 1 - e^^cos 2 ?jL + sin 2 tjL) - CON [e^^cos 277 L - sin 2 t?L) - ll 
— = _L i (4a) 

Kt ? e 2T?L sin 2r]L + CON (e 27?L cos 2r?L + l) 


where 


CON = A ( x , L' . ^Han <g - _ l and <p < 0 (4b) 

A(x, L, rj) + tan cp 

for the wall temperature lagging the fluid temperature 
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(4c) 


a(x, l, t?) = e 2 — 1 ^ 2L - x ) + cos m 

e 2rj(L-x) gin _ x ) + sin vpi 

where 



This solution is simplified if the wall temperature is measured at x = L, which for 
many applications is the easiest place to locate a sensor. The CON used in equation (4a) 
reduces to 

CON = tan (<p - r/L) at x = L (4e) 

The solutions to equation (1) may be used to determine the convective heat-transfer coef- 
ficient when thermal conductivity is an important factor. 

The computations made for the convective heat-transfer coefficient with the present 
solution (eqs. (4)) are time consuming. The values for h/Kry were calculated as a func- 
tion of 7 ] L and < p, when x = L, and are plotted in figure 2. Using x = L (temperature 
sensor located at the insulated surface) is reasonable because a temperature sensor is 
easily installed at this point. 

The ratio of the amplitude of the wall- to the fluid-temperature oscillation for x = L 
can be written by inspection of equation (A15). The ratio is 
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(a) r}l, 0.002 to 0.01. 

Figure 2. - Convective heat-transfer coefficient as function of frequency and phase shift with temper- 
ature measured at insulated surface. 
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Figure 2. - Continued. 
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Phase shift, (p, deg 
<c) tjL, 0.1 to 3.0. 
Figure 2. - Concluded. 
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and 


e G 


= tan 


1 



This ratio can be used to determine the magnitude of the fluid-temperature oscillation 
Tq that would yield a measurable wall -temperature oscillation at x = L. 

The determination of the effect of wall properties and plate thickness L on the cal- 
culation of h was aided by the expanding of equation (4a) to a series form at x = L. The 
resulting equation for hQ is 


h 


pCLq> 
-tan cp 


3 + tan% 2, 2 4 3, 3 

21 J] Lt + — Tj Li 

3 tan cp 3 


+ 5 + tan^ tp + 10 tan^ ^4^4 
5 tan^ <p 


(4g) 


In this equation, cp < 0. Convergence of this series must be checked. However, if 
rjL « 1 and reasonable values of cp (e.g. , -45°) are used, the series will converge. 

In this report, Anderson's solution is referred to as the slug solution because it 
equates the rate of change of the temperature of a mass or slug pL to the convective 
heat transfer from a fluid that flows over the slug. The slug solution is derived from the 
differential equation 


T 



+ T w " T G 


(5) 


and assumes that 

(1) The thermal conductivity K is infinite (i.e., no temperature gradient in the 

wall) 

(2) The convective heat-transfer coefficient h is constant 

(3) The density p and specific heat C of the wall are constant 

(4) One wall surface is insulated 

(5) The other wall surface is exposed to a fluid whose temperature is given by 

Tq = Tq + ATq sin it)6 (6) 

The solution to the differential equation (5) is 
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(7a) 


T = 


AT G 




where 


_ pCL 


h 


G 


(7b) 


is the time constant. The phase lag between the wall and the fluid is tan * wr = <p 
from which the convective heat-transfer coefficient can be written as 



-pCLo) 
tan cp s 


( 8 ) 


where w is < 0. 
Y s 


CRITERION FOR USE OF SLUG SOLUTION 

Comparison of the series form of the present solution (eq. (4g)) for 11 q with that of 
the slug solution (eq. (8)) shows that the coefficient of the series solution is simply hQ g 
If h r is substituted in the series solution and all terms having powers greater than 

Vj y S 

second power are neglected, the series solution can be written as 


h fi - = l - 1 + ton2 <P (r?L) 2 ... (9a 

h G,s 3tan< ? 

where <p < 0, and 77 L is assumed to be much smaller than 1. The value of <p can and 
should be approximately -45° (as discussed in the section Optimum Phase Angle). The 
curves in figure 2 can be used to determine the proper frequency that, for a given ma- 
terial, tjL, and h Q , will yield the value <p = -45°. The selected value of r?L is then 
used in equation (9a) along with (p = -45° to evaluate the ratio h^/hQ g . For these 
conditions, equation (9a) reduces to 



= 1 + -(hL) 2 
3 


= 1 + 2 
3 K 


G, s 


(9b) 
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where cp = -45° and 77 L « 1. The second term in equation (9b), an approximation to 

(h - h )/h, can be used to approximate the error in h if the slug solution, instead of the 
s 

present solution, is used to calculate h. The term may also be used as a first-order 
correction if the slug solution is used to calculate h. 

Values for the second term in equation (9b) were calculated for 347 stainless steel 
and copper and are compared in table I with the values of (h - h )/h x 100 percent. The 

o 

agreement is good, even though the absolute value of the difference is large in several 
cases. Care must be exercised in using equations (9a) and (9b) because only two terms 
of the series are used, and if the value of 77 L approaches 1, convergence is not ensured. 
In addition, values for cp that approach either 0° or 90° will greatly affect hQ/hQ . 

In addition to the aforementioned criteria to be used in the choice of solutions for 
calculating h, a number of calculations were made and are presented to illustrate the 
practical application of the solutions first given for h. 

APPLICATION OF EQUATIONS TO ENGINEERING MATERIALS 

The calculation of the relation between phase lag angle <p and forcing frequency f 
necessitates that assumptions be made for values of the convective heat-transfer coeffi- 
cient and the properties of the wall. The value of the convective heat-transfer coeffi- 
cient h^, assumed in this calculation is 0.001 Btu per square inch per second per °R 

w o 9 

(2.942x10 W/(m )(K)), unless otherwise noted. The wall properties used are those of 


347 stainless steel and are given as follows: 

Average temperature, T, °R; K 1000; 555.5 

Density, p, lb/in. 3 ; kg/m 3 0.286; 7.92X10 3 

Specific heat, C, Btu/(lb) (°R) ; J/(kg)(K) 0.128; 0. 536X10 3 

Conductivity, K, Btu/(in.)(sec)(°R); j/(m)(sec)(K) 0.000253; 18.92 


These values are approximately those of 347 stainless -steel materials used for simulated 
rocket-nozzle heat-transfer studies conducted in an air facility. The value of h^ is ap- 
proximately 5 to 10 times greater than those found by Anderson (ref. 1) in his wind-tunnel 
tests on a cone. 


Comparison of Phase Lags 

The calculated phase lags cp as a function of frequency f are shown in figure 3. As 
would be expected, the slug solution phase lag falls between the values calculated with 
equation (4a) (present solution) at x = 0 and at x = L. The decrease in ip at x = 0 for 
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Frequency, cps (Hz) 

Figure 3. - Phase lag as function of forced fluid-temperature frequency 
for 347 stainless steel. Fluid temperature, 1000° R (555. 5 K); convec- 
tive heat -transfer coefficient at fluid surface, 0.001 Btu per square 
inch per second per °R (2.942xl0 3 W/(m 2 KK)). 


f > 0.45 hertz can be explained by the reflected wave that counteracts the phase shift of 
the primary wave. 


Optimum Phase Angle 

When the phase lag is measured for use in either of the solutions, an optimum value 
exists. At this optimum value, a given error in phase lag will produce a minimum error 
in the heat -transfer coefficient. Figure 4 shows this optimum phase lag to be approxi- 
mately -45°, or exactly -45° if the slug solution is used. If equation (4a) is used with 
x = L = 0.040 inch (1.016X10 -3 m), the optimum <p increases slightly, which is the rea- 
son for choosing a frequency that will give a value for cp of approximately -45°. 

Phase Lag Differences 

The values for the phase lag calculated from both the slug solution and the present 
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0 20 40 60 80 100 

Phase lag, ip, deg 


Figure 4. - Change in convective heat-transfer coefficient for incremental change in phase 
lag as function of phase lag. (For present solution, temperature was measured at insu- 
lated surface. 


solution do not agree. The percent difference was calculated to determine the magnitude 
of this disagreement. The percentage is based on the present solution and is shown as a 
function of frequency in figure 5. For the case where the wall is 0.040 inch (1.016 

_ O 

X10” m) thick and the temperature measurement made is assumed to be at the insulated 
surface, the minimum difference is 6 percent. The frequency for this point is 0. 175 
hertz and the phase lag is -61. 5°. If it is possible to measure the temperature on the 
surface over which the fluid flows (x = 0), the difference can be reduced by using a lower 
frequency. For a phase lag of -45°, the difference at x = L is 6. 3 percent while at 
x = 0 it is 4 percent. 


14 




Figure 5. - Difference in phase lag due to location of 
temperature sensor as function of frequency for 
347 stainless steel. Fluid temperature, 1000° R 
(353.5 K); convective heat-transfer coefficient at 
fluid surface, 0. 001 Btu per square inch per 
second per °R (2.942x1Q3 W/(m 2 MK)) ; wall thick- 
ness, 0.040 inch (1.016x10'-* m). 


i 

i 



Frequency, f, cps, (Hz) 

Figure 6. - Difference in convective heat -transfer co- 
efficient due to location of wall temperature sensor 
as function of frequency for 347 stainless steel. Fluid 
temperature, 1000° R (555. 5 K); convective heat- 
transfer coefficient at fluid surface, 0.001 Btu per 
square inch per second per °R (2.942xl0 3 W/(mfy(K)) ; 
wall thickness, 0.040 inch (1.016xl0“ 3 m). 
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Heat -Transfer-Coefficient Differences 


Since a minimum value of the phase lag differences in the case where x = L was ob- 
served in figure 5, a similar minimum would be expected to exist when the heat-transfer- 
coefficient differences are calculated. Figure 6, however, shows that the lowest possible 
difference calculated for x = L is greater than 7 percent and occurs at a lower frequency 
than does the minimum phase lag difference shown in figure 5. At phase lags of -45°, the 
heat-transfer -coefficient differences are 9.6 percent at x = L and 7 percent at x = 0. 
The differences will increase with increased frequency or phase lag angle. 


Effect of Heat -Transfer Coefficient 


Up to this point in the calculation, the convective heat-transfer coefficient has been 



Figure 7. - Difference in convective heat-transfer co- 
efficient for 347 stainless steel. Fluid temperature. 
1000° R (355. 5 K) ; wall thickness, 0. 040 inch 
(1. 016xl0" 3 ml; temperature measured at insulated 
surface. 


O <3 0 

assumed to be 0.001 Btu per square inch per second per R (2.942x10 W/(m )(K)). 

Shown in figure 7 are the calculated heat-transfer -coefficient differences at x = L = 0.040 
inch (1.016x10"^ m) for h = 0.001 Btu per square inch per second per °R (2.942x10^ W/ 
(m^)(K)) and for h = 0.003 Btu per square inch per second per °R (8.826x10^ W/(m^)(K)). 
Figure 7 shows that, for low frequencies, the differences are greater for the higher h. 
When a phase lag angle of -45° is desired (to minimize the effect of errors in phase lag 
measurement), the differences are 9.6 percent for h = 0.001 Btu per square inch per 
second per °R (2.942x10^ W/(m^)(K)) and 23 percent for h = 0.003 Btu per square inch 
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per second per °R (8.826xl0 3 W/(m 2 )(K)). With a phase shift of -45° used as a criterion, 
it is concluded that increasing the heat-transfer coefficient will increase the heat- 
transfer-coefficient difference. 


Effects of Temperature 

The effect of wall temperature on the convective heat-transfer-coefficient difference 
is shown in figure 8. Because the temperature affects the wall properties, two tempera- 



Figure 8. - Difference in convective heat- 
transfer coefficient due to change in wall 
temperature as function of frequency for 
347 stainless steel. Convective heat- 
transfer coefficient, 0.001 Btu per square 
inch per second per °R (2.942x10^ 
W/(m%)) ; wall thickness, 0.040 inch 
(1.016X10" 3 m); temperature measured 
at insulated surface. 


tures were used, 1000° R (555. 5 K) and 2000° R (1101 K). The wall properties for 
347 stainless steel at 1000° R (555. 5 K) were given in the first part of this section and for 


2000° R (1101 K) are as follows: 

Density, p, lb/in. 3 ; kg/m 3 0.286; 7.92X10 3 

Specific heat, C, Btu/(lb)(°R); J/(kg)(K) 0.151; 0.632X10 3 

Conductivity, K, Btu/(in.)(sec)(°R); J/(m)(sec)(K) 0.000365; 27.29 


The values for the heat -transfer -coefficient differences calculated at x = L show that a 
2-percent decrease in the differences exists for a two-to-one change in wall temperature. 
Property changes due to temperature, in the case of 347 stainless steel, will not signifi- 

nonflu o Tf Aof f hrt « A+lfm Knof f r»n ncfor rlifforpnpoc 

Vb VllV CWilV 11V/1* l- V.J. 
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Effects of Wall Thickness 


Wall thickness can be expected to have an appreciable effect on the convective heat- 
transfer -coefficient difference, as shown in figure 9. For a frequency that will give a 
<p value of -45° and with the use of x = L, the following differences can be obtained: 
for L equal to 0. 010 inch (0. 254x10"^ m), 2.6 percent; for L equal to 0. 060 inch 
(1.524x10"^ m), 14.7 percent. For thin walls, the differences do not increase rapidly 
with frequency. The thicker walls cause the differences to increase at lower frequencies 
than those of the thin wall. Inspection of equation (9a) shows that the differences approach 
zero as L approaches zero. 



Figure 9. - Difference in convective heat -transfer coefficient due to change in wall thickness as 
function of frequency for 347 stainless steel. Fluid temperature, 1000° R (535. 5 K); convective 
heat-transfer coefficient, 0.001 Btu per square inch per second per °R <2.942xl0 3 W/(m 2 )(K)) ; 
temperature measured at insulated surface. 


Effects of Thermal Conductivity 

The effect of thermal conductivity on the convective heat-transfer-coefficient differ- 
ence is shown in table I. A comparison of copper and 347 stainless steel, having wall 

_ O O 

thicknesses of 0. 060 and 0.040 inch (1. 524x10 and 1.016x10 m), respectively, shows 

that the higher conductivity of copper decreases the heat-transfer-coefficient differences 
by 14. 2 percent. Equation (9a) shows that, as K approaches infinity, the differences 
approach zero. 
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COMPARISON OF SOLUTIONS 


The slug solution (ref. 1) , which neglects the effect of thermal conductivity and 
temperature-measurement location, may be used in place of the more complicated solu- 
tion (eqs. (4)) provided that the system is designed properly. Equation (4g) may be used 
to estimate the error in h when the slug solution is used provided that tjL < 1. 0. If a 
maximum error of 6 percent is to be tolerated, r/L cannot exceed 0. 2 and <p must be 
approximately -45°. Improper design will result in large errors. For example, a wall 

O 

made of 347 stainless steel, 0.060 inch (1.524x10 m) thick, with the temperature mea- 
sured at the insulated face will give errors greater than 23 percent if phase lags exceed 
-45°. Inspection of equation (4g) shows that if the slug solution is used, thin walls are 
essential. Although low frequency improves the accuracy of the slug solution, it is well 
to keep in mind that, at least for x = L, the limit of h^/h^ g ^ 1. Also, the accuracy in 
the measurement of the phase lag angle becomes very poor as <p approaches zero (see 
fig. 4) . High thermal conductivity is desirable as are low density and specific heat. If 
the phase lag is -45°, an increase in the heat -transfer coefficient will increase the error 
when the slug solution is used (see fig. 7) . An increase in the wall temperature of 
347 stainless steel from 1000° to 2000° R (555. 5 to 1101 K) resulted in only a 2-percent 
change in the error (fig. 8) . 

Table I summarizes the results of the comparison made of equations (4a) to (4e) 
and (8). This table presents calculations made for phase lags of -45°. 


CONCLUDING REMARKS 

The convective heat-transfer coefficient h can be calculated for a fluid flowing over 
a surface with one insulated side if the fluid temperature is varied sinusoidally. The 
phase lag between the fluid and wall temperatures, along with the frequency of oscillation 
and wall material properties, can be used to calculate the convective heat-transfer coef- 
ficient h. Two solutions for h are available. Both require a phase lag of approximately 
-45° to minimize the error in h due to errors made in measuring the phase lag angle. 
Anderson’s slug solution (ref. 1) does not account for the wall thermal conductivity or the 
location of the measured wall temperature, which may result in an error greater than 
23 percent in h. 
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A general one -dimensional solution is given which accounts for a finite thermal con- 
ductivity and for the wall-temperature location. This solution is greatly simplified if the 
wall temperature is measured at the insulated surface. Neither solution is applicable 
when two- or three-dimensional heat transfer in the wall is important. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 25, 1968, 

122-29-07-03-22. 
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APPENDIX - DERIVATION OF HEAT-TRANSFER COEFFICIENT 
AS FUNCTION OF PHASE LAG 


Determination of Boundary Conditions 

The temperature response of a wall, which has one surface insulated at x = L and 
the other surface exposed to a fluid with a temperature that varies sinusoidally at x = 0, 
is calculated as follows: First, the boundary conditions are determined with the assump- 
tion that 


T G ~ AT G e 


-iu)6 


For x = 0, 


£ 


G fc 


-iwd 



— 8T(0, 9) 

dx 


For x = L 


9T(0, 9) _ q 
9x 


The governing differential equation is 


£TM r l 9T(x, 6) 
dx 2 ' or d6 


For the solution to the differential equation, assume a product solution 

T = X(x) • F(0) 


(Al) 


(A 2) 


Then 


— = F(0) • X'(x) 
9x 


d T = Ffai 

— \~ f 

0X 


Y"(v\ 
— \—/ 
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and 


— = x(x) • F'(e) 

de 


Substituting these expressions in the differential equation (Al) gives 


F(e) • x”(x) = - x(x) • F(e) 

a 


or 


EMBLEM (A3) 

X(x) a m) 

Since either side of this equation is independent of the other variable, assume that each 

2 2 

side must be equal to a constant, X . Now, X can be equal to zero, greater than zero, 

2 

or less than zero. Then, setting either side of equation (A3) equal to X gives 

m = x 2 (a 4) 

X(x) 


and 

I ^ T ( g ) = x 2 
01 F(e) 

m = x 2 « 

F(e) 

Then 

In F = X 2 ad 

F = Ce X 01 6 (A 5) 

If X 2 = 0, F = 1 and the wall temperature T will not be a function of time. This result 
cannot be the case physically; therefore, the solution for X 2 = 0 is rejected. The 

p p 

choice between X < 0 or X > 0 is made by attempting to solve the equation by using 
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2 2 9 

-X and then +X , one of which will lead to a solution. With the use of -X , equation (A 4) 
becomes 


X"(x) + X 2 • X(x) = 0 


From Wiley (ref. 3, p. 88), the solution for this equation is 


X = Cje 1 ^ + C 2 e“ iAx 


(A 6) 


Equation (A 5) then becomes 


F 


= Ce 


-X 2 a-0 


(A7) 


o 

This equation is periodic when X is imaginary. The solution requires that the exponent 
be of the form c o9 . Therefore, X 2 is set equal to ico/a, and 


X = 



X = 



X = ±(1 + i) 'l/-^ 

2a 


Substituting for X and X 2 in equations (A 6) and (A 7) gives 


x = c e ±i(l+i) fiw/2a)x + Q e =Fi(l+i) ^(oo/2a)x 
1 2 


and 


F = Ce -1<00 

Substituting these solutions into equation (A2) gives a solution to the differential equa- 
tion (A3) which is periodic. 

Set 


V 
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The solution thus becomes 


T = + Be ±7?x e -i ^ 0:t ^ ( A8 ) 

where the constants A and B are CCj and CC 2 , respectively, and can be combined 
with no change in the solution. The choice of sign used in equation (A8) does not matter 
because the constants are arbitrary. The top sign is used in the following derivation. 

Applying the boundary conditions determines the constants A and B, but first some 
preliminary work is necessary. The boundary conditions require the use of dT/dx. 
Taking the derivative of equation (A 8) with respect to x gives 

H = ulAe-^e-^-^d - 1) + Be^ x e- i(w0+T ? x) (l - i)l (A9) 

0x L J 

Using boundary condition 2 and setting equation (A9) equal to 0 at x = L give 

A _ - e ^ e -iM+nx) (1 _ t) 

B e -r ? x e -i(ce0-r ? x) (i _ J} 

A _ e 2?]x-i2?7x 
B 

A = e 2rjL(l-i) (A10) 

B 

Using boundary condition 1 

aT G e" lwe - Ae- 1 " 8 - Be" 1 " 9 = ^[Ae^O - 1) - Be’^d - lj] 

h G L 


and dividing by 


e -iwe 


give 


AT g - A - B = ^2L[(A - B)(i - 1)] (All) 

h G 

Note that, if the fluid temperature is assumed to be T,-, = AT^e +ia;0 , the constants A 

ici)Q 

and B become functions of 6 because e cannot be eliminated by division (see eq. 
(All)). Solving equation (A10) for A and substituting in equation (All) give 
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AT^ - Be 27,L ( 1_i ) - B = - ~ K7 ? . 

Ci «, 

h G 


; e 2 ^<i-i) . B j ( . _ j 


Collecting terms and solving for B give 


B = 


AT, 


1 - e 2r 7 L ( 1-i )"J(i - 1 )*1 + [i e 2vUl-i) 
L J h~ L J 


Rearranging this equation gives 


B = 


AT, 


-Ki ? H _ e 2??L(l- ’M* 

^G 


Changing to the polar form gives 


B = . 


AT, 


-Kt, 


jl _ e 2r?L (cos 2r)L - i sin 2?}L)J + [l + e 2 ^(cos 2tjL - i sin 2 ? 7 L)j + i 52 jl - e 2?,Ij (cos 2 t;L - i sin 2?7L)j 

ho 


Collecting the terms in the denominator on i gives 


AT r 

B = 

(i - 5lA + /K2 + A e 2 ^ L cos 2r,L - 52 e 2r > L sin 2 t,L + i 

V h G / \ h G / h G 


-*n - l) e 2 ^ L sin 2?7L + — _ e 2r ' L 5H cos 2r?L 


G 


The following trigonometric substitution can be made in the previous equation: 



Kj 

he 
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where 


E = 


lt2 £U 2 (52 

h G \ h G/ 


Then 


AT G 

B = — 


1 - + Ee 27?L (cos £q cos 2tjL - sin e G sin 2r)L) + 

h G 


ISO. + e 2? ? L E(-cos e G sin 2rjL - sin e G cos 2?}L) 
h G 


This equation simplifies to 


B = 



1 - + e 2T?L E cos(e G + 2^L) + 

h G 


i ISH - e 2r?L E sin(e r + 2r]L) 
_ h G j 


The complex numbers must be in the numerator so that the phase shift accounts for the 
resistance of the boundary layer h G> This requirement will become apparent. To put 
the complex numbers in the numerator, divide the denominator into the numerator in the 
previous equation for B by multiplying each number by the conjugate of the denomina- 
tor. The following equation results: 


B = 


AT G e 


-U 


(A 12) 


n 


1 - ^ + e 2r?L E cos(e G + 2r]L) 
h G 


K^? _ sin(e G + 2-tjL ) 


L h G 


From equation (A 10), A is 


A = 


AT G e 2 ^ L e-^ +2 ^ 


(A 13) 


1 _ M + e 27]L E cos(e G + 2?}L) 


Kr] _ e 2rj L E sin(e r + 2pL) 


li?G 
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where 


£2 - e 27?L E sin(e G + 2r?L) 
h G 

£ = arc tan 

1-^2 + e 27?L E cos(e Q + 2 t 7L) 
h G 


From equations (A8), (A 12), and (A 13), the solution can be written as 
AT G ^( 2L -x) e -i(w0-?? x +l+ 2 77L) + e 77x e -i(co0+r?x+|)j 


T = 


(A 14) 


fl 


1 - 52 + e 27?L E cos(e G + 2?]L) 
hr- 


52 - e 2 ^ L E sin(e G + 2pL) 

hri 


If the real part of the driving temperature, that is, T G = AT G cos u>0, is selected, the 
imaginary part of the solution can be dropped. Equation (A 14) then reduces to 


T = 


ATj^le 77 ^ 2 ^ - ^ cos(a >9 - rjx + ^ + 2?jL) + e^ x cos(co0 + t?x + oj 


n 


(A 15a) 


1 - 52 + e 2r?L E cos(e G + 2?]L) 
h-, 


52 _ e 2?7lj E sin(e G + 2 17 L) 


where, as defined before but restated here, 


and 



(A 15b) 


52 - sin(e G + 2r?L) 

h r 

£ = arc tan — — (A 15c) 

1 - 52 + e 2??L E cos(e G + 2 t?L) 
h G 


o 

Equation (A 15a) is the required solution. The values for X >0 can be ruled out if the 
same process using X z is followed. Boundary condition 1 will then yield a solution for 
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2 

the constants that are functions of time, which cannot be; therefore, X > 0 is rejected. 

It should be stated that X 2 > 0 will work if T Q = AT Q e lw9 . The solution for this 
case will be the same as for equation (A15a). 


Determination of Phase Lag 

The phase lag is determined with a value found for cod such that T = 0. This con- 
dition will occur when the wall-temperature vector in the complex plane reaches -it/ 2 
radians. Therefore, the phase lag is 


cp = — - cod (A 16) 

2 

This quantity is less than zero for the wall temperature lagging the fluid temperature. 
Setting T = 0 in equation (A 15) gives 

e ??(2L-x) _ cos cod_ cos(tjx + Q - sin q>9 sin (tjx + Q 

e ?jx cos cod cos (-px + £ + 2tjL) - sin cod sin (-rjx + £, + 2tjL) 


Then 

e 2 ?? (L-x)j^ os (_ r?x + 4 + 2r;L) - tan cod sin (-?jx + £ + 2tjL)] = -cos(tjx + 0 + tan c od sin (tjx + |) 
from which, 


tan cod = e 2?? ( L ~ x ) cos(-tjx +J_+ 2tjL) + cospjx + j) 
6 (2 tj(L-x) + £ + 2tjL) + sin(Tjx + |) 


and 


6 2tj(L-x) C0S H: + 7 j(2L - x)| + cos(| + tjx) 
cod = arc tan — “ — 

e 2rj(L-x) gin ^ + tj( 2L - x)] + sin(£ + tjx) 


(A 17a) 


The heat -transfer coefficient can then be determined from the phase lag. Solving 
equation (A 17) for £ and then for hQ gives 
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(jo9 = arc tan e^ L ~ x ) [cos g cos r?(2L - x) - sin g sin tj(2L - x)] + cos g cos 77 X - sin | sin tjx 
e 2^(L-x) j" s j n ^ cog - x) + cos | sin ?j(2L - x)] + sin £ cos tjx + cos £ sin 77 X 


e 2r?(L-x) tan ??(2L _ x) ] + _cosj7x _ . ten 4 __sinjx 

0,0 = arc tan cosj^L^x) c os tj ( 2L - x ) 

e 2 ?7( L -x)n n 5 + tan tj(2L - x)l + tan | ^osjx + _§!“JE 

L j /OT /OT 


cos tj( 2L - x) cos tj(2L - x) 


0)0 = arc tan 


-tan £ 

e 2r?(L-x) - x) • 

sin 7jx . „2 tj(L-x) . cos rjx 

cos ?7(2L - x) cos 77(2L - x) 

tan | 

' 2 t)(L-x) COS TJX 

" 1 c 2t?(L_x) tan 7i(2L - x) + sin ^ 

cos 7j(2L - x) 

T C Will 1 jy£j JLJ aj t 

I cos ?}(2L - x) 


(A 17b) 


From equation (A 16), 


o >9 = -- <p 

2 


where cp is less than zero for the wall temperature lagging the fluid temperature, and 

tan o)0 = tan - (pj 
tan o)0 = cot cp 

The foregoing expression is used in equation (A 17a) to write 


cot cp<tan £ 


J^2t7(L-x) 


COS TjX 


L 


cos ?7(2L - x)J 


+ e 2 ”< L - x > tan „(2L - x) + s » n ”* 

COS 77(2L - x) 


+ tan 4 


|g2f7(L-x) 


L 


tan tj( 2L - x) + 


sin T7X 


cos t?(2L - x)J 


- e 


2tj(L-x) cos tjx _ q 

cos tj( 2L - x) 


Solving for tan £ gives 
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(A 18) 


e 2T 7 (L- X )r 1 _ ten v{2L _ x)cQt I + cos jjk _ sin rjx cot^ 

L J ~ /Ox /OT 


tan | = 


cos r/( 2 L - x) cos t?( 2 L - x) 


>( L ' x) + C0S ^ l oot <p + e 2 ^ L - x ) tan r,( 2 L - x) + 

cos r/(2L - x)J cos 77(2 


■ S CON 


?7X 

r^( 2 L - x) 


This expression is set equal to CON for simplicity and the definition of £ is used to 
write 


- e 2T?L E sin(e r , + 277 L) 


CON = tan | s 




1 - — + e 2 ^^E cos(€q + 2rjL) 
ho 


The solution for h^K rj is as follows: Keep in mind the trigonometric substitution for E 
and €q and write 


K V . e 2r?L 


_ 




|- 


^cos2„L + (l + ^ 
_ h G V h G) 

| sin 277L 

= CON 4 

l.K3 t e 2 " L 
[ h G 

[1 + !S?\ cos 2r?L - -V sin 2??L 
_\ h G/ h G 

► 

> 


|i . e 2?7L (cos 277 L + sin 2 tjL) - CON [e 2T?L (cos 2rjL - sin 2r?L) - ljj 

= e 27?L sin 2rjL + CON(e 2r?L cos 2r?L + 1) 

Solving for IIq/Ktj give 

h Q 1 - e 2T?L (cos 2r]L + sin 2r,L) - CON[e 27?L (cos 2?]L - sin 2r]L) - l] 


—± (A 19) 


Kr? 


e 2??L sin 2r?L + CON(e 2t?L cos 2r;L + 1) 


Simplifying equation (A 18) by first multiplying the numerator and denominator by 
tan cp, and collecting terms gives 
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CON =_ 


LWL-x) tan „ (2l . x) + sinyx — 1 + r e 2ij(L-x) + cos „x ^ 

L cos r^(2L - x)J L cos 7?(2L - x)J 


+ r e 277 (L-x) + cos rpt 

cos ?]( 2 L - x). 


- + e 2 ^](L-x) ten _ x ) + 


sin 77X 


cos 77(2 L, - x)J 


tan cp 


Defining 


A(x,L, 77) 


6 2t7(L-x) + cos 77X 

cos 77(2 L - x) 

e 2? 7( L -x) tan t?( 2L - x) + 

cos 7](2L - x) 


A(x, L, 77) = 


e 2r ?( k~ x ) cos t^( 2L - x) + cos rjx 
6 2 tj(L-x) g ^ n - x) + sin 77X 


Then 


CON = L ’ ^) tan £ JL.. 1 (A 20) 

A(x, L, 77) + tan cp 

Equations (A19) and (A20) are the solutions presented in the text. For x = L, 

A(x, L.,77) = cot 77L 


CON _ tan <p - tan 77L 
1 + tan cp tan 77L 

CON = tan(<p - 77L) 


(A 21) 
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